function [i_line] = cal(t)
    i_line = -U0_C_current/L_current/wd*exp(-alpha*t)*sin(wd*t);
end
